Topology, phase transitions and the spherical model 
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The Topological Hypothesis states that phase transitions should be related to changes in the 
topology of configuration space. The necessity of such changes has already been demonstrated. We 
characterize exactly the topology of the configuration space of the short range Berlin-Kac spherical 
model, for spins lying in hypercubic lattices of dimension d. We find a continuum of changes in the 
topology and also a finite number of discontinuities in some topological functions. We show however 
that these discontinuities do not coincide with the phase transitions which happen for d > 3, and 
conversely, that no topological discontinuity can be associated to them. This is the first short 
range, confining potential for which the existence of special topological changes are shown not to be 
sufficient to infer the occurrence of a phase transition. 
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Phase transitions (PTs) remain one of the most intrigu- 
ing and interesting phenomena in physics. Mathemati- 
cally, a PT is signaled by the loss of analyticity of some 
thermodynamic function pj in the thermodynamic limit. 

Recently, a new characterization of PTs has been pro- 
posed, that conjectures that "at their deepest level PTs 
of a system are due to a change of the topology of suit- 
able submanifolds in its configuration space" J3|. This is 
known as the Topological Hypothesis (TH) |3J]. In this 
new method one studies the topology of the configuration 
space r of the potential energy V(x) of a system with N 
degrees of freedom, determining the changes that take 
place in the manifolds M v = {x e T : V(x)/iV < v} 
as the parameter v is increased. A topological transition 
(TT) is said to take place at c if M c _ e and M c+e are not 
homeomorphic. The idea is that somehow TTs may be 
related to PTs. 

The necessity of TTs at a phase transition point has 
been demonstrated for short ranged, confining models 
In the XY model Q TTs are found both in the mean 
field (MF) and unidimensional short range (SR) versions, 
whereas a PT is present only in the MF case. This led 
to a refinement of the TH: only sufficiently "strong" TTs 
would be able to induce a PT. It was found that in the 
MF version a macroscopic change of the Euler charac- 
teristic happens at exactly the same point v c where a 
PT appears. Several other models seem to be in agree- 
ment with this behavior But recently it was proved 
for a nonconfining potential that no topological criterion 
seems to be sufficient to induce a PT We show below 
that the same happens for the spherical model, which is 
a confining, short ranged potential. 

In the spherical model it has been found that there 
is a direct correlation between the TT and the PT, in 
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its mean field version Q. Interestingly, in the case of 
nonvanishing external field there is no PT, but the con- 
figuration space displays a TT at energies that cannot be 
thermodynamically reached. 

In this letter we study the original Berlin-Kac spheri- 
cal model [8] for spins placed on a d-dimensional lattice, 
interacting with their first neighbors. Using tools from 
topology theory we were able, for the case of vanishing 
field, to determine its topology exactly (up to homology). 
We show that the PT occurring for d > 3 cannot be re- 
lated to any discontinuity in the homology of the mani- 
folds at v c . For nonvanishing field we cannot character- 
ize the topology completely for all v, but show that a 
very abrupt change in the topology happens that does 
not have a corresponding PT. At variance with the MF 
version, the value of v at which this topological change 
occurs is thermodynamically accessible. 

The spherical model is defined by a set of N spins e, ly- 
ing on a d dimensional hypercubic lattice and interacting 
through the potential V = —\ J2<ij> Jij e i e j " H J2i e i 
where Jij = J gives the strength of the interaction be- 
tween nearest-neighbor spins i and j, and H is an ex- 
ternal field. The spin variables are real and constrained 
to lie on the sphere S N ~ X (i. e. J2i e i = N). Periodic 
boundary conditions are imposed on the lattice. 

In H it is shown that, at zero field, a continuous PT 
appears at a critical temperature T c (d) for d > 3, which 
is a strictly increasing function of d (see Table 1). On 
the other hand, no PT is possible in an external field. 

As in previous works, the thermodynamic function we 
use to relate the statistical mechanical and topological 
approaches is the average potential energy per particle 
(v). Although the specific details of (v) depend on d, 
some features are common to all hypercubic lattices: 
(v) — ► for T — > oo and (v) — > — d (its lower bound) 
when T — ► 0. 

In the topological approach one looks for changes in the 
topology of M v as v is increased. A topological change 
happens at a certain value vt if the manifolds M VT ~ e and 
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3.9573 


-1.0216 
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6.4537 


-0.7728 
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8.6468 


-0.6759 
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10.7411 


-0.6283 



TABLE I: Critical temperatures T c and mean potential ener- 
gies per particle (v c ) for hypercubic lattices in d dimensions. 
Values obtained from analytical expressions in Q. 



M VT+e are not homeomorphic Q for arbitrarily small e. 
To make a connection with statistical mechanics Casetti 
et al. 0] proposed the nontrivial ansatz that, at the phase 
transition, vt can be identified with (v c ), the thermo- 
dynamical average critical potential energy per particle. 
To study the topology of the configuration space of the 
spherical model it is most convenient to write the po- 
tential using the coordinates xi that diagonalize the in- 
teraction matrix through an orthogonal transformation: 
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where we set J = 1, and A^ (i = 1, • • • , N) are the eigen- 
values of the interaction matrix, ordered from largest to 
smallest. We define the sets Cj, j = 0, • • • , N, where 
N + 1 is the number of distinct eigenvalues. Cj is the 
set containing the indices of the eigenvalues that have 
the (j + l)-th largest value. Therefore, \Cj\ gives the de- 
generacy associated to the (j -I- l)-th largest eigenvalue. 
The Frobenius-Perron theorem ensures that the largest 
eigenvalue is not degenerated, i. e. Co = {1}. 

The critical points of this potential on the sphere 
T = S^- 1 = {x G R N : J2f=i x2 j = N ] are found us- 
ing Lagrange multipliers. Along with the spherical con- 
straint, the critical point equations are: 



:ri(2^ + Ai) + VNH = 

Xi {2n + Xi) = i = 2,- 



(2) 
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where fi is the Lagrangian multiplier that results from 
enforcing the spherical constraint. From these equations 
and Eq. N + 1 critical values of v are obtained, de- 
noted Vk = — A;/2, with I £ Cfe, and k = (),-•• ,N (or- 
dered from smallest to largest). Notice that the degener- 
acy of the eigenvalues causes that the corresponding crit- 
ical points be in fact critical submanifolds. This implies 
that in the directions tangent to the critical submani- 
folds the Hessian vanishes, which in turn implies that the 
potential is not a proper Morse function. Nevertheless, 
using Bott's extension of Morse theory the Euler char- 
acteristic can be found exactly. However, profiting from 
the symmetries of the spherical model we took a more 
direct route to study its topology. As we show below, 



for vanishing external field it is possible to characterize 
completely the topology of the M v , by explicitly giving 
the values of all the Betti numbers of the manifolds. 

For H = the critical manifolds £„ . , j = 0, • • • , N, 
are given by T, v . — {x S T : X«ec 3 x i = N } ( see ec L- 
l|2jl). These are (hyper) spheres whose dimension is given 
by the degeneracy of the corresponding eigenvalues. To 
understand the nature of the topological change that hap- 
pens at the critical values of v it is necessary to know the 
topology of the M v for v between two critical values. We 
show below that in the interval (vj,Vj + i) all the mani- 
folds M v are homotopy equivalent to S D_1 , where D is 
the number of eigenvalues larger than —2v. In fact we 
prove that S 13-1 is a deformation retract of M v , which 
in turn implies their homotopy equivalence 0|. 

A submanifold S C M is a deformation retract of a 
manifold M if there exists a series of maps f v : M — > M 
with v e [0, 1], such that f° = I, f l {M) = S and f\ s = 
I for all v. The map when considered as / : M x [0, 1] — > 
M must be continuous. Let us take v £ {vj^Vj+i). The 
deformation retract that takes the manifold M v onto its 
submanifold S 13-1 = {x 6 M v : x 1 — N} is given 

byxM = (/ 1 "(x) ) ... ,/ft(x)) with ~ 
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for i < D 



for i > D 
(3) 

This map can easily be shown to be continuous at all 
points x € M v . The properties for v = and v = 1 are 
evidently fulfilled. It is also easy to see that the retraction 
f v does not map any points outside M v , since the image 
points always lie on the sphere S^ -1 , and their potential 
energy does not exceed v. 

Homotopy equivalence implies that the Betti numbers 
of the M v with v € (vj, w^+i) are the same of S 13-1 : 
h(M v ) = 1 for i = and i = D - 1 and h{M v ) = 
otherwise. Thus at each Vj a topological transition oc- 
curs that changes the topology of the phase space from 
one homotopy equivalent to S 0- '^' -1 to one homotopy 
equivalent to S D_1 . In terms of the Betti numbers, each 
transition changes two of them from to 1 and from 1 to 
0. Thus, at variance with other models, the magnitude 
of the Betti numbers is not a useful quantity in order 
to characterize the TT. It is better to look at changes 
in D — 1, the highest index of the Betti number that 
changes at each transition. Furthermore, as we have 
shown that the manifolds M v are homotopy equivalent 
to (hyper) spheres, the information about their dimension 
D — 1 completely characterizes their topology. Thus D 
is the relevant quantity to be studied. As shown above, 
the increase of D at each TT is given by the degeneracy 
of the corresponding eigenvalue. If the degeneracy \Cj \ is 
o(N), given that D = £- = o |Ci|, in the N -> 00 limit D 
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is equivalent to the order of the critical manifolds, which 
is defined as the number of negative eigenvalues of the 
Hessian (J^iZo when restricted to the submanifold 
normal to the jth critical submanifold. This generalizes 
to degenerate manifolds the definition of order of a saddle 
point. 

For the spherical model it can be shown that the spec- 
trum of eigenvalues is continuous in the infinite N limit. 
Thus, the set of N + 1 critical energies will be dense in 
[— d, d], the interval of allowed potential energies. Con- 
sequently the model has a continuum of TTs. In this 
limit, and considering that D is O(N), it is convenient 
to introduce a continuous and normalized version of D, 
d(v) — D/N, and also a degeneracy density, c(v). They 
are related by c(v) = dd g^ ■ In the following we search for 
singularities in these functions or their derivatives which 
could point to particularly strong TTs. 

The spectrum of the adjacency matrix is given by ||: 
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PIG. 1: Degeneracy density c(v) and relative index d(v) of 
the critical manifolds, for a vanishing external field. 



A p = 2^cos(2 m /7V 1 / d ), Pi = <),■■■ ,N 1 ' d -l (4) 



In the N — > oo limit, the degeneracy density is: 



c(v) = (2tt) 
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where A(w) = 2 J2i=i c °s(wi). It can be shown 0] that 
the integral converges uniformly for all values of d and 
therefore c(v) is a continuous function (see Fig. QJ. The 
derivatives with respect to v can be obtained by perform- 
ing the derivative inside the integral. But, as this is only 
valid if the resulting integral converges, this procedure 
allows us to obtain only the first [(d — 1)/2J derivatives. 
All these derivatives are continuous [ll| except for the 
last, which is discontinuous only at the following points: 
at odd values of v if d is odd, at even values of v/2 if 
d/2 is odd and at odd values of v/2 if d/2 is even. But 
these values are clearly different from the ones at which 
a PT takes place, for all values of d (see Table 1). The 
manifolds M v display TTs at the points (v c ) where PTs 
occur, since there is a continuum of TTs. These TTs, 
however, are not particularly abrupt. This non coinci- 
dence between the levels where a special TT (vt) and a 
PT ((v c )) take place has also been observed in the <j) 
mean field model [l^ . 

The only possibility left to look for a relationship be- 
tween TT's (in the sense of a discontinuity of some func- 
tion of the topology) and PT's would be in the higher 
derivatives of c(v), which cannot be studied by inter- 
changing the integral and derivative operations. This 
possibility seems to us rather unreasonable, because it 



would imply not only that the derivative where disconti- 
nuities are to be looked for depends on the dimension of 
the lattice, but also that those discontinuities present in 
lower order derivatives should be disregarded. 

We have thus shown that discontinuities in the deriva- 
tives of c(v) are not sufficient to induce the PT present 
in the model. Furthermore, we show in the following 
that, even though in the case of a nonzero external field 
there appear discontinuities in the function c(v) itself no 
connection between such TTs and PTs can exist, simply 
because the model does not display any PTs at all. 

With H it is not so easy to find the homotopy type 
of the submanifolds M v , because of the breaking of the 
symmetry introduced by the field term in the Hamilto- 
nian. Nevertheless, using Morse theory it is at least pos- 
sible to establish the homotopy type of the submanifolds 
up to above the second smallest critical energy, where an 
abrupt topological change is shown to take place. 

According to Morse theory, if there is one nondegen- 
erate critical point at c € (a,b), the manifold M c+e is 
homeomorphic to M c _ e U eft, where is a k-ce\\. In other 
words, at the critical point, a fc-cell (i.e. a fc-dimensional 
open disk) is attached to the manifold, where k is the 
index of the critical point, defined as the number of neg- 
ative eigenvalues of the Hessian at that point. 

From the critical point equations (2J we obtain that the 
smallest critical energy is v+ = — (Xi/2+H), and the next 
is V- = — (Ai/2 — H). The Hessian of the potential on 
the sphere at the critical points x± = (±vJV, 0, • • • ,0) 
is a diagonal matrix with V u = X\ ± H — \ (i > 1). 
Therefore, at these two points the Hessian is not singular, 
which implies that x± are nondegenerate critical points. 
But Ai is the largest eigenvalue, therefore for v+ all the 
eigenvalues Vu are positive. This was to be expected 
because this is the absolute minimum of the potential. 
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Topologically this means that for v < v_ M v is homotopy 
equivalent to a disk. For v = v_ the index of the critical 
point depends on the field: if A 2 < Aj — H, w_ is a 
minimum. Thus, denoting the next critical value by v 2 , 
M v for v 6 («_, 1)2) is homotopy equivalent to the union 
of two disjoint disks on the sphere. However, for large 
values of N the topological scenario is diferent. Since in 
this limit the spectrum of the adjacency matrix becomes 
dense, a certain number k of its eigenvalues will fall into 
the interval This number becomes the order of 

the critical point at w_ , and gives the dimension of the 
fe-cell that is attached to the disk. The manifold M v for 
v e (v- , V2) is therefore homotopy equivalent to a sphere 
of k dimensions. For large N, k becomes proportional 
to N. In the interval the manifolds M v have 

the homotopy type of a point. At V- an abrupt change 
in the topology takes place, and the M v have now the 
homotopy type of a sphere with a macroscopic number 
of dimensions (see the jump of c(v) in Fig. 2). 
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for all values of v. Nevertheless we have found exactly 
the topological change that takes place at f_, and have 
shown that it is macroscopic. It may come as a surprise 
that this very abrupt change does not have a PT associ- 
ated to it. 

We have shown that the manifolds of the configura- 
tion space of the short range spherical model display a 
continuum of topological transitions. Hence the neces- 
sity condition implied by the theorem in @j is trivially 
met. Also, strong discontinuities have been found either 
in a function of the topology or in its derivatives. Al- 
though these discontinuities represent abrupt changes in 
the topology we have shown that they are not associated 
to PTs. Conversely, at the points where PTs take place 
no abrupt changes are observed in the topology. These 
are the first results on a short range confining potential 
to challenge the sufficiency of a topological mechanism in 
the origin of a phase transition. 

One intriguing question that arises is whether the topo- 
logically abrupt changes that we have found can have 
some influence on the dynamics of the model. 

This work was partly supported by CNPq (Brazil). S. 
R-G. acknowledges support from the Centro Latinoamer- 
icano de Flsica. 
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FIG. 2: Degeneracy density c(v) and relative index d(v) of 
the critical manifolds, for a finite field. A delta function is 
assumed at the discontinuity point in the graphs of the left 
column. 

For higher values of V, the critical values are given by 
vj = —\j/2 + H 2 /2(Xi — Xj), but only for j such that 
Ai — Xj > H . Notice that, at variance with the case 
of vanishing H, there is a threshold energy below which 
the critical values have been suppressed. The critical 
submanifolds occuring at each Vj are again hyperspheres 
whose dimension is given by the degeneracy of the cor- 
responding eigenvalue Xj. For all critical values we have 
calculated the order of the critical manifolds, d(v), as well 
as the relative degeneracy, c(v) , for a few values of d (see 
Fig. 2). The main difference with the results for H = 
is that now the connection between the order and the 
topology of the different manifolds is less obvious, and 
we have not been able to identify the homotopy types 



t Electronic address: srisau@cab.cnea. gov. ar 
* Electronic address: anacarol@if.ufrgs.br 
^ Electronic address: stariolo@if.ufrgs.br Research Asso- 
ciate of the Abdus Salam International Centre for Theo- 
retical Physics, Trieste, Italy 

[1] C. N. Yang and T. D. Lee, Phys. Rev. 87, 404 (1952); T. 
D. Lee and C. N. Yang, Phys. Rev. 87, 410 (1952). 

[2] L. Casetti, M. Pettini, E. G. D. Cohen, Journal of Sta- 
tistical Physics 111, 1091 (2003). 

[3] L. Casetti, M. Pettini, E. G. D. Cohen, Physics Reports 
337, 237 (2000). 

[4] R. Franzosi and M. Pettini, Phys. Rev. Lett. 92, 
060601 (2004); R. Franzosi, M. Pettini and L. Spinelli, 
math-ph/0305032 

[5] R. Franzosi, M. Pettini and L. Spinelli, Phys. Rev. Lett. 
84, 2774 (2000); L. Angelani, L. Casetti, M. Pettini, G. 
Ruocco and F. Zamponi, Europhys. Lett. 62, 775 (2003); 
P. Grinza and A. Mossa, Phys. Rev. Lett. 92, 158102 
(2004) 

[6] M. Kastner, Phys. Rev. Lett. 93, 150601 (2004). 
[7] A. C. Ribeiro-Teixeira and D. A. Stariolo, Phys. Rev. 
E 70, 16113 (2004). See however an alternative per- 
spective on the problem by M. Kastner, preprint 
cond-mat/0412199 
[8] T. H. Berlin and M. Kac, Phys. Rev 86, 821 (1952). 
[9] A. Hatcher, Algebraic Topology (Cambridge University 
Press, Cambridge, 2002). 
[10] R. Bott, Bull. Amer. Math. Soc. 7, 331 (1982). 
[11] S. Risau Gusman, A. C. R. Teixeira and D. A. Stariolo, 

to be published. 
[12] A. Andronico, L. Angelani, G. Ruocco and F. Zamponi, 
Phys. Rev. E 70, 041101 (2004); D. A. Garanin, R. 



Schilling and A. Scala, Phys. Rev. E 70, 036125 (2004). 



